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INTRODUCTION  AND  REVIEW  OP  LITERATURE 

The  free- convection  effect  in  vertical  pipe  flow  has  recently 
become  quite  significant  In  a  number  of  engineering  applications. 
The  modifier  "vertical''  here  is  employed  to  mean  that  the  longi- 
tudinal dimension  of  the  pipe  is  parallel  to  the  direction  of  the 
generating  body-force,  such  as  gravity  or  centrifugal  force.  The 
two  important  fields  where  the  free- convection  phenomenon  plays 
a  dominating  role  nowadays  are  haat- exchanger  problems  in  nuclear 
reactors  and  blade-cooling  problems  in  high  speed  turbines. 

The  field  of  pure  free- convection  and  combined  free-  and 
forced-convection  pipe  flow  are,  as  to  be  expected,  extremely 
barren  compared  to  that  of  pure  forced-convection.  Contrbutions 
of  a  theoretical  nature  are  confined  to  those  for  par all el- plate 
channels  and  circular  pipes.  Ostrach  analyzed  pure  free-convection, 
constant  heat- genera ting  laminar  flow  with  constant  wall  temper- 
atures (10);  combined  free-  and  forced-convection,  constant  heat- 
generating  laminar  flow  in  channels  with  linearly  varying  wall 
temperatures  (11)  and  unstable  flow  in  channels  with  wall  temper- 
atures linearly  decreasing  in  the  vertical  upward  direction  (12 
and  13).  Ostrach  also  pointed  out  that  energy  dissipation  can  be 
neglected  only  when  a  parameter  (the  product  of  the  thermal  coef- 
ficient of  volumetric  expansion,  the  body  force  per  unit  mass,  the 
hydraulic  diameter  of  the  pipe  cross  section  and  the  reciprocal  of 
the  specific  heat  capacity  at  constant  pressure)  is  small.  The 
effect  of  energy  dissipation  has  been  investigated  by  successive 
approximations  in  all  Ostrach' s  works  except  one  (13).  Lietzke  (*+) 


-• 


analyzed  pure  free- convection  laminar  flow  without  heat  generation 
in  channels  with  both  ends  closed,  using  parameters  different  from 
those  employed  by  Ostrach.  In  Lietzko's  analysis,  heat  fluxes 
through  the  walls  are  prescribed  and  energy  dissipation  was  neg- 
lected entirely, 

Lighthill  (5)  investigated  pure  free- convection  in  circular 
pipes  with  one  or  both  ends  closed.  Both  laminar  and  turbulent 
flows  were  considered  with  energy  dissipation  neglected,  Lighthill' s 
analysis  is  the  most  inspiring  among  the  only  three  (2,  5  and  1*f) 
that  deal  with  the  effects  of  the  presence  of  the  pipe  ends.  A 
method  of  obtaining  solutions  satisfying  the  integrated  forms  of 
the  governing  equations,  which  distinguishes  Lighthill1 s  work  from 
the  others  reviewed  here,  has  been  employed  by  him,  Hammitt(2),  and 
Ostrach  and  Thornton  (1*+)  extended  further  the  laminar  part  of 
Lighthill* s  problem  by  the  same  integral  method, 

Hallman  (1)  treated  combined  free-  and  forced-convection,  heat- 
generating  laminar  flow  in  circular  open-end  pipes  with  linearly 
varying  wall  temperature  in  a  way  similar  to  Ostrach* s.  However, 
in  Kallman's  treatment,  energy  dissipation  was  neglected  and  differ- 
ent parameters  were  used.  Furthermore,  as  quoted  in  Ostrach  (11), 
approximate  solutions  of  co  bined  free-  and  forced- convection  flow 
in  short  pipes  and  channels  were  obtained  by  Kartinelli  and  Boelter, 
Works  by  Ostroumovj  Hanilton,  Poppendiek,  Redmond  and  Palmer; 
Hamilton,  Poppendiek  and  Palmerj  Woodrow;  iiurgatroydj  Wordsworthj 
and  Timo  were  also  quoted  by  Ostrach  (11  and  12)  and  Hallman  (1). 


They  all  covered  the  same  aspects  as  Ostrach  (10,  11  and  12), 
Lietzke  (k)   and  Mailman  (1),  Through  private  conjmini  cation,  Dr. 
E.  R.  G.  Eckert  informed  the  author  that  Lietzke  has  extended 
his  work  on  parallel-plate  channels  to  the  case  of  circular 
annuli  but  the  result  has  not  been  made  accessible. 

Izuni  (3)  solved  unsteady  free  connection  inside  vertical 
pipes  by  Laplace  transform.  This  work  is  made  known  only  through 
Ostrach' s  review. 

The  unrestricted  way  of  attacking  the  problem  of  the  thesis 
would  be  to  solve  itwfltimiilji  the  continuity,  the  -Javier- 
Stokes  ,  and  the  energy  equations  for  a  compressible  fluid  with 
body  force  and  energy  dis  ipation  terms;  aided  by  the  equation  of 
state  and  the  laws  of  variation  of  the  specific  heat  capacities, 
thermal  conductivity  and  dynamic  viscosity.  But  all  the  existing 
solutions  are  for  se:;;i-incompressible  fluids  with  constant  specific 
heat  capacities,  thermal  conductivity  and  dynamic  viscosity.  The 
fluid  is  considered  to  be  generally  incompressible,  but  a  buoyancy 
term  is  introduced  by  taking  into  consideration  the  variation  of 
density  with  respect  to  temperature.  This  trend  is  also  followed 
in  the  present  investigation.  A  general  theory  of  the  present  pro- 
blem, applied  to  pipes  of  any  cros  section,  is  developed  first. 
This  is  essentially  a  generalization  of  Ostrach' s  work  on  parallel- 
plate  channels.  Gases  where  the  axial  component  of  the  ter.perature 
gradient  vanishes  are  set  apart  from  those  where  it  does  not.  In 
all  the  existing  solutions,  algebraic  elimination  of  temperature 


from  the  energy  and  the  Navier-Stokes  equations  is  always  employ- 
ed. In  the  present  investigation,  however,  algebraic  elimination 
of  velocity  is  also  used,  as  sometimes  this  is  more  convenient, 
especially  when  applying  thin-plate  analogy.  To  the  best  of  the 
author's  knowledge,  this  thin-plate  analogy,  discovered  in  the 
course  of  the  present  investigation,  has  not  yet  been  discussed 
in  the  literature.  This  analogy  is  not  essentially  the  same  as 
that  in  pure  forced-convection  problems  (.arco  and  Han,  6),  since, 
generally,  elastically  supported  plates  serve  as  models  here. 
Other  useful  analogies  are  also  used. 

Following  the  general  theory  thus  established,  the  problem  is 
solved  for  a  few  specific  cross  sections.  Flow  inside  circular 
annuli  is  treated  fully.  Solutions  for  rectagular  pipes  and  pipes 
of  ring- sector  cross  sections  are  obtained  by  using  the  finite 
Fourier  sine  transform.  This  and  some  other  mathematical  means  are 
introduced  as  fomal  tools  only.  Mathematical  rigour  is  not  attempted. 
The  above-raentioned  specific  solutions  are  the  main  contributions 
of  the  present  investigation  and,  so  far,  are  not  available  in  the 
literature.  Investigations  regarding  parallel-plate  channels  and 
circular  pipes  have  been  done  rnther  extensively  by  others,  and 
hence,  are  not  duplicated  here.  Some  special  cases  for  triangular 
cross  sections  are  solved  by  applying  thin-plate  analogy.  Possibi- 
lity of  solution  by  other  analogies  or  techniques  Is  pointed  out. 
Ho  attempt  is  made  to  follow  Ostrach  in  investigating  the  effects 
of  energy  dissipation.  The  phenomenon  of  unstable  flow  is  not  explor- 
ed. It  Is  also  to  be  noted  that  only  the  problem  of  prescribed  wall 


temperature  is  taken  up.  The  possibility  of  solving  the  eases  of 
elliptical  pipes  and  annuli  (composed  of  two  confocal  ellipses) 
in  terms  of  Kathieu's  functions  and  other  possible  directions  in 
continuing  the  present  investigation  are  pointed  out  in  the  SUGGES- 
TIONS, 

GENERAL  THBORY 

General  Assumptions 

General  assumptions  mentioned  explicitly  or  implicitly  in  the 
INTRODUCTION  AND  REVIEW  OF  LITERATURE  are  summarized  here  as  follows t 

1.  The  pipe  wall  is  cylindrical. 

2.  The  elements  of  the  pipe  wall  are  parallel  to  the  generat- 
ing body  force. 

3.  The  domain  investigated  is  far  away  from  the  ends  of  the 

pipe. 

k*   The  flow  is  steady  and  laminar. 

5.  The  fluid  is  seml-incompressiblej  i.e.,  f(see  Appendix  I 
for  nomenclature)  is  treated  as  constant  except  that  the  relations 

P-PCT) 

1  dp 

and  |38--- 

'     p  dT 

are  used  when  taking  the  buoyancy  effect  of  density  variation  into 

consideration.  In  practice,  suitable  values  of  p,  a,  and/or  ^3  may 

be  used. 

6.  c  ,  c  ,  k  and  u  of  the  fluid  are  constants. 

7.  Energy  dissipation  is  negligible. 
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8.  dT/dz=con3tmt=4,  i.e.,   T=Az+T*(x,  >x2). 

9.  The  temperature  distribution  on  the  vail  is  known. 

Coordinate  System  and  Domain  of  Variables 

A  general  Cylindrical  coordinate  system  is  used  (Fig.  1).  A 
Garti  ian  coordinate  axis  z   is  chosen  parallel  to  the  elements  of 
the  pipe  vail  and  negative  in  the  direction  of  fz.  The  mid-point 

of  the  length  of  that  portion  of  the  pipe  where  the  flow  is  well 
developed  is  chosen  as  the  origin  of  z.  In  any  plane  perpendicular 
to  the  z-axls,  a  two-dimer:oional  coordinate  system  (x..,x2)  is  es- 
tablished. Surfaces 

x<j=constant 
and  X2=constant 

are  cylindrical  x/ith  elements  parallel  to  the  z-axis.  For  convenience, 
any  plane  perpendicular  to  the  z-axis  will  be  referred  to  as  the 
two-di-ensional  region,  as  differentiated  fro;:  the  three-dimen- 
sional domain  consisting  of  the  planes 

and  the  cylindrical  pipe-wall  surface  y  which  is  the  boundary  of 
the  two-dimensional  region  8« 

General  Equations 

Under  the  assumptions  r.entioned  above,  the  continuity  equation 
becomes  simply 

dw 


—  ^=0 
dz 


which  i  plies  that 
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v  ■  w(x-j,x2),  only. 
Another  useful  form  of  the  continuity  equation  is 

pllvdCT  ^constant  mass  through-flow  =g  (1) 

S 
where  the  mass  through- flow  g  is  a  known  parameter  of  the  problem. 
The  first  two  Navier-Stokes  equations  give 


(grad  p)  =0 


x1 


and  (grad  p)  =0 

*2 

which  imply  that 

p=p(z),  only. 

The  third  Navier-Stokes  equation  assumes  the  form 

jUV^  -  —  -pfz  =0.  (2) 

dz 

The  buoyancy  effect  can  now  be  brought  into  equation  (2)  by  adding 

and  subtracting  a  reference  density  termj 

/v\ .  (p_Pri)fz .  | ,  prA  (3) 

where  Pr  is  p   evaluated  along  an  element  C,  off.  Now  assume 

and  use  the  definition 

1  dp 
'     (°  dT 
Three  different  cases  must  be  considered  separately  here  (Ostrach, 
10): 

1.  If  temperature  differences  in  the  whole  region  are  small, 


. 


p  and  p  are  essentially  constants  and  suitable  average  values  of 
them  can  be  used  such  that 

or  -<p-pr>)=p(3(T-Tro). 

2.  If  temperature  differences  are  large  but  the  fluid  is  a 
liquid  obeying  the  relation 

{>p>  =constant, 

dlJ  PmPr 

then.  —  =  -Pfl  =constant=  — U*. 

dT    '  r  T-fr, 

or  -<P-?r->»pP<*-*r.>- 

(If  f|3  is  not  exactly  a  constant,  suitable  average  value  of  it 

could  be  used.) 

3.  If  temperature  differences  are  large  but  the  fluid  is  a 
gas  obeying  the  lav 

-  J- 
BT 

with  p-pr  (so  as  to  be  consistent  with  the  assumption  that  P  is 
a  function  of  T  only),  then, 


I     la  rr,  tvti 


r°  3Tr 


|MrJ 


BT^T 


=  JL(T-Tr  ),  or  -S-u(T-Tr  ) 


10 


But 


I3 

1   dp       1 
?  dT       T  ' 

h; 

■  i/*r.  * 

or 

Here  suitable  average  values  of  p  or  p  can  be  used  so  that 

For  all  these  cases  equation  (3)  becones 

provided  that  suitable  average  values  of  P,  (5  and/or  pp  are  select- 
ed for  each  specific  case.  0  in  equation  (k)   stands  for 

T-Tfa  =  (Az+T*)-(Az+T*o  )=T*.T*c)  , 
and  is  a  function  of  x^  and  x2  only.  The  left-hand  side  of  equation 
(h)   is  a  function  of  x1  and  x2  only  while  the  right-hand  side  is 
a  function  of  z  only.  Therefore,  both  sides  of  equation  (k)   must 
equal  to  a  constant,  say  II,  and 

juv2v  +  ppfze  =11.  (5) 

For  an  incompressible  fluid,  the  energy  equation  gives,  neg- 
lecting energy  dissipation, 

pCyAw  *  Q+  kv2!  (6) 

[noting  that  du  =(du/dT)vdT  +(du/dv)Tdv  ■  c^dT,  since  dv=o]  or, 
since  Tro  is  a  function  of  z  only, 

f>cvAw  -   Q  +  kV  e.  (7) 

The  key  equation  of  the  problem  can  be  obtained  by  ellminat- 
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ing  •  from  equations  (5)   and  (7)* 

v\  +  f  f  ?°v  A  w  =  C£j6  q.  (8) 

^k  yUk 

When  boundary  conditions  are  specified  equation  (8)  can  be  solved 

for  v;  e  can  then  be  obtained  from  equation  (5).  The  parameter  II 

involved  in  the  solution  can  be  expressed  in  terms  of  the  known 

mass  through-flow  g  using  equation  (1).  In  this  way,  it  is  not 

necessary  to  explicitly  define  the  type  of  flow.  Whether  a  flow  Is 

pure  free-convection  or  combined  free-  and  forced-convection  is 

implicitly  embodied  in  the  numerical  value  of  g. 

The  problem  can  also  be  solved  by  eliminating  w  from  equations 

(5)  and  (7)  i 

,    o2Af  c  A    Pc  A     \t2Q 

v\  ♦  r  ft  g  v  e  .  Lx.  ii  .  _.         (9) 

yis.  jJk  k 

After  e  is  obtained  from  equation  (9),  equation  (7)  gives  w.  The 
parameter  II  can  again  be  expressed  in  terms  of  g  through  equation 
(1). 

If  it  is  proved  to  be  more  convenient,  equatios  (8)  and  (9) 
may  also  be  paired  together  to  give  w  and  e  without  resort  to  equa- 
tions (5)  and  (7). 

Boundary  Conditions 

There  are  two  kinds  of  principal  boundary  conditions.  One  of 
them  is  the  non-slip  condition  of  a  viscous  fluid  at  the  wall J 

w  =0  on  T  (10) 

or,  using  equation  (7), 

V  ©  = on  f.  (11) 

k 
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The  other  is  the  thermal  condition  at  the  wall  (which,  in 
the  present  investigation,  is  the  temperature  distribution) J 

•  =  »p  on  p 
=0   at  fo 
Using  equation  (5),  equation  (12)  can  be  written  as 

.2    1 


(12) 


7^  --(II  -f|Sfz«r)  onr 


II 


(13) 


Beside  these  two  kinds  of  principal  boundary  conditions  there 

is  a  number  of  auxiliary  ones  which  might  be  needed  at  times,  e.g., 

that  w  and  o  are  finite;  or  that 

-oo<g<e>o  when  both  ends  are  open, 

g=0     when  one  or  both  ends  are  closed. 

When  one  or  both  ends  are  closed  there  must  be  pure  free- convection 

in  the  pipe. 

Special  Cases  with  A=0 

When  A=0  equations  (7)  and  (8)  may  be  written  as,  respectively 

kv2©  +  Q  =0  (1*0 

v*v  *  llLs  q.  05) 

Equations  (1),  (5),  (10),  (12)  and  (13)  still  hold.  Equation  (9) 
becomes  the  same  as  equation  (1*0.  Equation  (11)  holds  in  the 
whole  region  S,  instead  of  on  J"  only,  and  thus  becomes  the  same 
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as  equation  (1*0. 

Here  the  problem  may  be  solved  by  determining  •  from  equation 
(1*0  and  then  w  from  equation  (5)  J  II  is  expressed  in  terms  of  g 
by  equation  (1).  It  may  also  be  solved  by  finding  v;  from  equation 
(15)  and  then  e  from  equation  (5),  or  w  and  e  may  be  solved  inde- 
pendently from  equations  (15)  and  (I*!-). 

Nondimensionali  zation 

In  order  to  nondimen  ionalize  all  the  equations  mentioned 
above,  the  hydraulic  diameter  d  of  3   is  used  as  the  characteris- 

/c  ra 

tic  length.  For  cases  with  Aj*Of   /  A  i —  and  Ad  are  used  as  the 

p     Pr 

characteristic  velocity  and  temperature,  respectively.  For  those 


with  A=0,  /-" — "-  and  Tv     are  used.  In  this  manner,  dimension- 
P      Pr        lo  ' 

less  equations,  grouped  together  for  quick  reference,  are  obtain- 
ed as  follows | 
A*0 

GA=jfAdr  (16) 

S 

k 

V  «A+RaA«A»«A  (17) 

^A=0  on  T  (18) 
T2^A=FA-rAr  on  f  ■) 

=FA           at  J  (19) 

1  ▼V  VFA  (20) 


1»f 


TjfV  on  r  | 

T2^  -  §LA  on  r 


RaAHl*^ 


W* 


T« 


^=0  on  f 

▼%fFA  -  %• on  r  | 

*FA  at  ro- 

TA=^r  on  f 
=0       at  To 

*2V  -  ** on  r 


A=& 

G=\Ldr 

s 

tJ=0  on  r 
/  ▼^0=P  -tr    on  T 
«F  at  Tc 


) 


(21) 

(22) 

(23) 
(2*0 
(17) 
(18) 

(19) 

(21) 

(22) 
(23) 

(25) 

(26) 
(27) 

(28) 
(29) 
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^T2?-*  ^  =0 

(30) 

tsTTp  on  r  "* 
|     =0     at  roJ 

■ 

(3D 

T2co+  f  «F 

(29) 

^co»0  on  r 

(27) 

^^♦■0;  =  ^ 

(26) 

cob  o  on  T 

(27) 

{ 

▼2co«F  -tr  on  T  \ 
*F             at   To  J 

X*\  on  p  > 

^    =0     at  r0J. 

(28) 
(30) 
(3D 

In  rectangular  Cartesian  coordinate  system, 

d» 

«•  d£d7 

(32) 

T2 

d2 

d2 

(33) 

* 

a2 

• 

tm 

In  polar  coordinate  system, 

d* 

*  RdRd0 

(35) 

T2 

16         6 
R  dR       dR 

1 

R2 

d2 
"dP 

(36) 

* 

1   1 

=( fl 

R  dR 

6 

I  —  + 

to 

1 

— r)2 
602 

(37) 
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Thin-Plate  Analogy 

An  inspection  of  equation  (21)  shows  that  it  is  analogous  to 
the  equation  governing  small  deflection  of  a  thin  plate  resting 
on  an  elastic  foundation  with  TA,  RaA  and  RaAFA-r  iA  correspond- 
ing to  the  deflection,  the  modulus  of  the  foundation  divided  by 
the  flexural  rigidity  of  the  plate,  and  the  lateral  load  divided 
by  flexural  rigidity,  respectively  (Timoshenko,  17,  f  H?)«  If  the 
boundary  conditions  happen  to  be  analogous,  too,  an  exact  similarity 
exists.  :A:  ilar  comments  apply  to  equation  (17). 

Equation  (26)  shows  an  analogy  to  the  equation  governing  small 
deflection  of  a  thin  plate  supported  at  the  edges,  without  elastic 
foundation  lying  beneath  it  (Timoshenko,  17,  ^21). 

In  applying  the  thin-plate  analogy  and  borrowing  known  solu- 
tions from  that  field,  so  e  generality  must  be  sacrificed,  especial- 
ly when  velocity  is  compared  with  deflection.  In  the  following,  this 
analogy  is  introduced  only  when  the  more  general  cases  are  not 
attempted.  Other  analogies  are  also  mentioned  at  the  end  of  SPECIFIC 
SOLUTIONS. 

Heat  Transfer 

The  Nusselt  number  used  in  this  investigation  Is  defined  as 

hd   qd     it  i      &T-  I 

Hu  =  —  0  =  -  — d /  9° m r%  (38) 

k   e°k   dN  /      dW 


■ 
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or 4  /  T;  (38) 

d  W   A 

where  ©°is  a  certain  reference  temperature  difference  and  Tr"  or 

t"A  is  its  dimensionless  counterpart.  The  selection  of  «°  is  a  matter 

of  convenience  only,  as  '/ill  be  shown  later, 

SPECIFIC  SOLUTIONS 

Circular  Annuli 

In  this  specific  case,  T  is  composed  of  two  concentric  circles 
with  diameters  I\-AD,.  and  Dj,  Since  the  hydraulic  diameter  is  Dx-^jt 

the  outer  and  inner  dlmenslonless  radii  are 

1 

1    *"1  I  (39) 

A-1  , 
respectively. 

Further  assumptions  are  made  that  ft  and  T^  are  two  constants 

and  that  Q  is  independent  of  0,  Then  the  problem  becomes  axLally 
symmetric,  I.e.,  d/d0  =0  and  the  following  equations  may  be  ob- 
tained, Tj  being  chosen  as  T^  t 

Ml 

SR7V 
^ARdR  (*K)) 

hi 

1  d    d 


,» — «— : 

H  dR   dR 


( R  — )^wA  +RaAcoA=  &A  (Vl) 
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°°A=0  when  H=Rj  or  % 

1  d         d 
( r  — )w    =p    when  R=fi_ 

flda       dB     A      A  * 


"Va  when  RSi^ 


1   d         d 
( R  — ) 

RdH       dR 


A+     A=FA 


t 
H 


G=2nl     WRdR 


Ri 


(^2) 


(V3) 


(W 


W) 


1   d         d 

/     R  dR       dR 
1^=0  when  ftlftg  -s 

=(1-^)F  when  R=RA  J 

1   d         d 

( R  —  )co+  "tr    =F 

R  dR       dR 

o»*0  when  Mb  or  R^. 

The  solutions  of  equations  (*K))   to  (¥0,  when 

1   d         d 
( *  — }d\  =0, 

R  dR       dR       A 

are  (Appendix  II) 
1  r 

<°A  =  -[^^^(s^^^beioCsARj+A^keroCs^O^V^o^A11)]*  VSA  CJW 


(*f6) 

(^7) 

0*8) 
(**9) 


(50) 


S?  r 


A  =  FA  +     7  [A1bGlo(SAa)-^bero(sAR)+^keio(sAR)-Vero(sAH)] 


1     1   d         d     * 

+  -t< ■  — )A. 

s  7  R  dR      dR     * 
A 


(52) 


„ 

^^{^^•(^)^^<^^)]^[1^<v^^•A^ 

+A3 

RXk^^sAHA)-Veio(sAHI)]A[HAkero(sAHA>-HIkero(sA^ 

27U*. 
SA  "I 

(53) 

where 

bero(sARI)        beio(sARI)       ^tWV       keio(sARI) 

A  « 

bero(sARA)        beio(£AR  }        keroUAHA)        kelo(sAHA> 
be^Cs^)     -fcer^Rj)       kei^s^)     -ker^s^) 

bei0(sARA)     -ber0(sAR  )       kei0(sARA)     -ker0(sARA) 

\M 

VSA       ^o^aV        kero(sA%)        keio(sARI) 

1 
SA 

VSA       belo(sARA)       kero(sAHA)       kelo(sAHA> 

FA       -^o^aV        ^o^aV     -kero<   AHI> 

Vl       "  bero(sARA>     k<*o<sARA>     -kero(sARA> 

bero(sABI)        VSA       kero(sARI)       kelo(sARI) 

(55) 

1 

v-7 

SA 

bero(sAR^)       VSA       kero(sAR>)       keio(3ARA) 
belo(sARI)           FA        kelo(sARI)     -ke^Cs^lj) 
beio(sARA)        eAFA        kei0(sARA)     -ker^s^) 

bero(sARI)       beio(sARI)       VSA       keio(A) 

(56) 

A3=.~ 
SA 

ber0(sARA)       beiQ(sARp       iA/sj[      kei0(sARA> 
bei^s^)     -ber^s^)           PA       -ker^s^) 
beiQ(sARA)     -berQ(sARA)        ^       -ker^s^) 

ber0(sARj)       beiQCs^j)       ker^s^)     %/sA 

(57) 

1 
A  = 

bero(sARx)       beio(sARA)       ker^s^)      i/s* 

beloCaAHI)     -bero(sARI}       kelo(3ARI)         FA 
^i0(sAR^     -bero(sARA)       keio(   ARA>      %FA 

(58) 
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Theoretically,  FA  in  equations  (51)  to  (52)  must  be  expressed 
in  terms  of  GA  with  the  aid  of  equation  (53).  But  it  is  always  easi- 
er in  practice  to  work  backwards,  i.e.,  to  assume  a  value  for  P. 
and  clculate  the  corresponding  G.. 

For  the  calculation  of  Mu,  Ad  is  chosen  as  «  ,  so  that 

dTA     i  dr 

and  Nu  T  *• S  =  -  — 2. 

i»A         dZ/        IdR 


A1  bel0<  sAR)-A2bero(  sAR)+^3keio(  SAH)  Akero(  sAR)j 

1     d     1   d         d      .  ■) 

♦  -r( R~ )^Ai(59) 

s*  dH  R  dH       da     AJ 

Both  Nuj  and  Nu^  can  be  calculated  from  equation  (59)  using  Rj  and 

R^,  respectively.  The  negative  sign  is  for  the  inner  wall}  the 
positive  for  the  outer. 

In  case  equation  (50)  does  not  hold,  the  particular  integral 
of  equation  (1*1)  depends  on  the  form  of  Q  and  is  not  simply  <i/s*J. 
In  practice,  fortunately,  Q  Is  either  zero  or  a  non-zero  constant, 
so  that  equations  (51)  to  (53)  always  apply.  For  simple  forms  of 
Q  the  particular  integral  may  be  obtained  by  standard  methods. 

It  is  of  interest  mathematically,  however,  to  note  that  equa- 
tion (1+1)  can  be  solved  without  any  restriction  on  Q  by  a  special 
finite  Hankel  transform  developed  by  Sneddon  (16).  The  result  is 
(Appendix  III) 
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*■'?■£*- *k 


A"  4 "         ww  ?!.♦■* 

(60) 

where  the  sun  Is  extended  over  all  the  positive  roots  of  the  trans- 
cendental equation  of  ^ 

Jo(Sa%)Go<^A)-Jo^mHA>Go<rniRl)^ 
With  equation  (60),  the  whole  problem  is  considered  to  have  been 

solved  mathematically.  But  the  solution  in  this  form  is,  of  course, 
not  easily  manageable  when  numerical  results  are  needed. 

For  special  cases  with  A=0,  equation  C+6)  yields  straight- 
forwardly the  solution 


-Jtf 


7r=  _l  -(^RdR)dR.  (61) 

a 

The  two  integration  constants  are  to  be  evaluated  by  using  equation 

0*7). 

In  case  Q  is  a  constant, 

t«  -  +CilnR  +C0# 

Using  equation  (*+7),  C-j  and  G2  are  determined  and 

A   >+1 
(1-£)P  +  — ' 

^«  .  -(R2-R2)+  ln(R/RT).        (62) 

h  T  In  A  X 

Equation  (51*)  gives,  generally, 

<*B  J  ~  ]J  [P  +3  -KpRdR)dRj  RdR/dR.  (63) 

Equation  0*9)  can  be  used  to  determine  the  integration  constants 
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involved. 

When  Q  is  a  constant, 


4  if 

co=  — R  +•£■ 


*        L     1+ 


i  0  16  A-1     0 

p-   -R^   +   R^ 


(1-6)  F+ 


k  X+1 
16  A-1 


InA 


^lnCR/Rj) 


InA 
+C1lnR  +  G2. 
C1  and  C2  can  easily  be  found  by  using  equation  (^9).  The  final 
result  is 


cox  JL<A.R^)+-$- 

61+    l 


i  Am 
(1-OP  ♦ 

F  -  -Rf  +  (R2-R?) 


(1-OP  + 

+  ^ (R2-R2)ln(R/RT) 

InA       A        * 


1  f  i   A*-1     1 
lnAL102»f  (A-lT   16 


3/  o 
F  -  -Rf  + 
If1 


£  A+1 
(1-OF  + 

16  ^-1 


InA 


yin—  (6*f) 

A-1J   Hj. 


The  egression  for  G  when  Q  is  constant  can  be  obtained  by  a 
straightforvaru  integration  of  equation  (6k).   The  result  is  too 
lengthy  to  quote. 

For  the  crlculation  of  m,  ©A  is  chosen  as  ©°,  so  that 


and 


NuI,\"  + 


1   d'fr 


When  Q  is  constant, 


Ku_    =    + 


H*     A(1-0F 


,A      J  (1-6)P  dR" 


i   A+1 
(1-€)F  +  JL 

i  16  A-1   1   | 


ln^ 


RJ  . 


(65) 
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Note  that  when  Hj-*-0,  a  circular  annulus  does  not  approach  a 
circular  pipe,  but  a  pipe  with  a  thin  wire  along  its  axis,  as  a 
limit. 

Rectangular  Pipes 

In  this  case,  F  Is  composed  of  four  straight  lines:  x=0f  x=a, 
7=0  and  y=b=o<a<.  Since  the  hydraulic  diameter  is  2ab/(a+b),  the 
dimensionless  side-lengths  are 

!*i(^+1)  (66) 

b  3  i     (1+0i) 

Choosing  T(0,0)  as  Tpof  the  following  equations  are  obtained: 

m 

b  a 

GA=)  \  °idjdy  (67) 

0  0 
d2        d2 

f  {7f +  of  )2wa  +  HSaWa  ^a  (68) 

coA=0  when  ^»0,   P=a,  »J=0  or  »j=b  (69) 

d2         d2 


■r  TA3(^)  when^a 

~fa"  ta2^^  when  7*° 

SFA~    ^^  When  f**' 


(70) 


*2         a? 


\  ,d  6 

1    '   -I 


V*J3H^A«*1  (71) 


£& 


0  0 


§d7 


d2    b2 

T«  t1(,7)  when?  =0 

■  't*M  when  <jr«a 

/     «  T2CO  when  ^=0 

=  "\(F)  when  t?«B 

62         d2 
(—-  +  — -)co+T=F 


(72) 


(73) 


(71*) 


\ 


d^2      d72 


(75) 


'ca)=:0  when  F=0,   Fsa,  viaso  or  w=b.  (76) 

Equation  (68),  with  the  boundary  conditions  (69)   and  (70),   can 
he  solved  by  double  finite  Fourier  sine  transform  (Sneddon,   16). 
The  result  is  (Appendix  IV) 

*        r  ir  i  ib       na 


00     00 


*?         r  _ir  1  m"       "«• 

af.p    i-(.i)»    i-(-On  (— +  — ) 

A     AL  -"-  J    nn        mb 


,     m2        n2   _      I 


mTYf         n?r>7 

•sin-— ^sin—^ 
a  b 


?[V«>nal*  ?&■">"* 


n2 

-o)  +sT 


a        b 


•sin — >sin — [(77) 
a  b 


The  second  double  series  in  equation  (77)  nay  be  simplified  to 
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CO 


CO 


-16F 


V  /  ,    ei2     n2  _ 

a^     b'-        A 


1     n2     n2 
( + ) 

nn  I2     b2  nm?       n7v? 

sin — -sin— — -• 
a  b 


e»«      »o 


m^      117*7 
sin — *sin— r~ 


k  m         n  ,p     k 
a         d 


a 


frO 


e»o 


-16F 


1     m2       n2 

b  ida|      nr^ 


A/  /  m2       „2  Sil 

m=1,3,5  n=1,3,5  7>  (~  v  ~>  +SA 


a 


BlTr 


VW)"V 


2  2 

nr       n^ 


sin—sin— 
a  b 


r[v(-1>"T*tl 


0"-       n^ 


m>F      nr**? 
-sin-— sin-^r*  • 
a  b 


ac 


*         A 


(78) 


Usually,   T^T^T^T^TCOjO),   i.e.,    \^^K\f^^  *  tha* 
equation  (78)   becomes 
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PO       c*a 


I 


a2  r    A 


~  in-  ■':-ln  _  ' 
a  b 


oo 


oo 


-16F, 


1     in        n 

una         b 


■ 


2  2 

»*,m        n  *2L> 


BiTtF       n-rv*? 

da — -sin-Tr*  .  (79) 


fT^5  ^7375  ^  ♦  r^+sj 

a         b 


If  in  addition,   Q  is  constant, 

1 


OO 


oo 


co  =16 


rai 


n2       n2 


ZA  _p  (—  +  — ) 

Lt^       A^2       ^'j 


m> 


2  2 

et       n^ 


=17375  fefT575  ££  +  3»*»»; 


-sin «  • 

b 


(80) 


Equation  (80)  in  turn  gives 

1 


00 


00 


W16 


icn 


-.     m2       n2  m2       n2    "> 

Va2       B2  A  a2       b2    - 


ePT73T5  n=i,j,^ 


2 

a*1       n' 


■sin 


rr 


**(_  +  — )2+slf 


00 


x> 


=16 


BSD 


fer^^rr^    ^%*^**4 


\    m2       n2       s;    i 

si  ( +  — -)+  -Af 

-  A  a2       b2      n2  AJ       m-rr?      n-M7 

-sin— -  sin L 

a  b 


m2       n2 


. sin—^ 

b 


(81) 


(81«) 


00 


00 


16F*\  \  1         ia-rt?      nrry 

(since  FA=  — -TV  /  — sirt — ^sin-^  .) 

*    I -  i -  m        a         b 

m=1,3,5  n=1,3,5 
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OC 


and  G.= 


6^fab 


re 


_         1      rL  m2      n2 

m2n2  Ltx2  A  a2       b2  . 

2         2 

m  n 


(82) 


/ , , ....  /, u,  ni        n    o    h. 

m-1,3,5  n=1,3,5      tT^  +  ~2>  *SA 


a^       W 

For  the  calculation  of  Nu,  Ad  is  chosen  as  e°    so  that  equations  (38) 
and  (81)  give 


CO 


CO 


2     2 
7,  '       m     it 


2  »2     ^ 

2 


an  L  A  a2  b2  A  ?2  >,2 


1  dtA 

5  bf1,3,5  n=1f3f?       ttV^  *  ~>  +SA 


a^  b< 


-si 


n— =H 
b 


CO 


c>o 


7T 


2      2 
m*  n^ 


m2  n2  j7 


2  drA 


mb  L  A  a2  b2  A  a2  B2 


m2       n2 


0=1,3,5  n=1,3,5       >fcj  +  =2)2+sA 


a 


a^       b' 

It  is  to  be  noted  that  equation  (81»)   can  not  be  differentiated 
term  by  term  since,   in  expanding  F&  into  its  double  sine  series, 
a  finite  Jump  from  -FA  to  Fa  Is  introduced  at  £=0,   J=a,   >^=0  and 
7=b. 

Similar  techniques  can  be  used  to  solve  equations  (73)  to 
(76).   The  results  are  as  follows: 

b 


k     \      \  a   L  1  JJ     b    L  ^  ^J       ETtf      n^r? 

~"  T?  (      I       5 5 sin sin— r~* 

ab  I /  „  nr       n^  a  b 


a         D 


(8*0 


#      m?rr 


60« 


ab  /        ■  j,  m 


^-<-i)mt3]+  2TjT^.(.i)^ 


m=1   n=1 


2  r,2 

a"       b2 


m7r5       rt7\w 
sin— ^sin— r-' 
a  b 


CO 


16 


m2      n2 


^~  • 


m=1,3,5  n=*1f3,5  anU.  +       ) 


a' 


When  T^^-Tn-T^^O  and  Q  is  a  constant, 

1 


oo 


oo 


16a 


ri7t 
rr-sin — =%! 


f        nTH7 

-sin— < 


rt 


■ 


=1,3,5  n=1,3,^  mn(-r  +^) 

a       \r 


n2     n*  a 


CO 


16 

7 


h         m2      n2 
00      ±  -F(—  >  —) 
tt2        a2 


b*1        m^j      nTry 
_ iia  _  tia  g "I 

m=1,3,5  n=1,3,5  on(~  +  ~)2 


ac 


b2 


oo 


oo 


m2       n2 


G=6»f 


aV 


L  .f(-~  ♦  — .> 

7T2         a2       b2         m?^      n7*7 

- —— sin— -^sin-r r*  • 

m2       n2  ^         a  b 


10=1,3,5  n=1,3,5  m2n2(—  +  -r)2 

a^       b^ 


(8if») 


(85) 


(86) 


(87) 


All  these  equations  are  also  derivable  fror    those  for  A^O  by  drop- 
ping out  the  subscript  A  and  making  sA  equal  zero.  For  the  calcula- 
tion of  Nu, 


oo 


a      b      165 

2       2?^ 


(.ir"-(-D 


n-1 


2         2 
,  m         n 

n=1,3,5  n=1,3,5  mn(— +  —3} 

a^       b 
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m+n 
. 22— - 22—  -1 

TFL  L n2       n2 

m=1,3,5  n=1,3,5  mn(-r»  +  ^) 

a        b 


is  chosen  asf°»  and 

CO  CO 


m+n 


«*1        _*S__£2_  — -1 

7t\  V  I "     /\  \  (-D 

Nu 


1»3       a/ / e2       n2  /  / / m2       n2 

m=1,3,^  r,=1,3,5  n(— ■  +  ~)     m=1t3,?  n=1,3,5  ninC--  ♦  ~) 

a2       b2  /  a2       b2 


(88) 


m+n 


i     °°       ,     *°  sin^      1,      °°      ■     °°  — -1 

"V**^ / m2       n2    // / m2      n2 

m=1f3t5  n=1,3,5  ra<— -  +  — )     m=1f3^  n=1,3,5  mnC^r  +  r~) 

a^       b*  /  a*1      b^ 

(89) 
Pipes  of  fling-Sector  Cross  Sections 

In  this  case,  p  is  composed  of  tvo  radial  lines,  0=0  and  0=§, 
and  two  circular  arcs  with  diameters  Dj   and  D^=ADi.   Since  the 
hydraulic  diaiaeter  is 

A2-1    : 

d=  ,  (90) 

A+1 

A-1+  | 

2 


the  dimensionless  radii  are 
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1      1 


(91) 


and 


V 


+£ 


<A+1)$   71-1 

Note  that  for  a  full  sector  Dj-0  and  A  =°%  while  Dj  remains  finite 
and  equals  to  D.  Hence 


d= 


tfAD, 


-L-H-— 

1  1 

H  «   1-  - 


1+if 


(92) 


tad 


lljSQ 


RA= 


(93) 


.1 


i   2 

When  $=2^  the  ring-  or  full  sector  becomes  an  annulus  or 

circle  plus  one  radial  line.  When  $  is  small  enough,  the  sector 
may  be  used  to  approximate  a  narrow  isosceles  trapezoid  or  tri- 
angle. 

Choosing  T(Ri,0)  as  Tr  ,  the  following  equations  are  obtained: 

Mi 


wAadRd0 


(9*0 
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1   a     d    1     a2  P 

(95) 

<^=0  when  R=Ri,   R=Rx ,  0=0  or  0=§ 

(96) 

< 

1     d        d       1       b2 

( r  _  +  _.  — _)u)  =F  -  T.,(0)  when  R=R  -> 

R  dR       dR       R2  d02     A     A       M                          1} 

=FA-  t^ (0)  when  R=RA 

=F  -  T»t^R)  wiien  0=° 

> 

(97) 

-F  -  tjx(^)  when  0=$  ' 

\  1     d        d      1       d2 

R  dR       dR  +  R2  d02  ^       A~  A 

(98) 

Ml 

f    I 

G=\    \        RdRd0 

(99) 

V% 

1     d         did2. 

( a  — -  +  -r  — *yt+sf=o 

/     R  dR       dR       R^  d0^ 

(100) 

T«srtI(0)  When  R=Rj  -\ 

< 

=  T*(0)  when  R=R^   I 
=  ^(R)  when  0=0     | 
=  T*(R)  when  0=§    ' 

(101) 

( 

1b         did2 

(-  --  r  --  +  —  -~)oHr«F 

R  dR       dR       R2  d02 

(102) 

u*=0  when  R-Rj,   R^R^  0=o  or  0=$  • 

(103) 

Equation  (95) »  with  the  boundary  conditions  (96)   and 

(97),   can 

be  solved  by  finite  Fourier  sine  transform.  The  result  is 

(Appendix 
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oo 


CO. 


n=1 


Aj 


n?r0 


A    5 


(1M 


where  ^A  is  the  solution  of 
2  d3<£.   1 


d  CO, 


,,,22  a2a. 
2n  -tt  a  co» 


2n2^  d 


^ 


dR 


R  dR- 


Rc 


§2  'dR2 


1     2n&7^  dal 
R- 


P  P  P  P 

r1  n  -Tr  n  7T 

+L7  ,2  <  X2  -'t)+Ha' 


^0   •, 


subject  to  the  boundary  conditions 
/CO^=0  when  R=Rj  or  R=R;\ 


3 

1  n?rf. 
R2  § 


§2 


°a=V  3  r|["-<-1>n]VV(-1>nW(105) 


S 


d2^.   1  dcq   q    r 

_A  + 4  =  L.FAl1-(-1)n 


dR' 


R  dR 


n7t 

I 

nor 


~T     when  ML 


-^'aI1-*-""]-^— »»««. 


For  special  cases  of  full  sectors,  Rj=0  and  t&£~0.   The  boundary 
conditions  for  equation  (105)  become 

Aj 

'U>k=0   when  R^O  or  R=R^ 


d26o.  1  dw. 


<? 


-~&  + **  —  PA  1-(-1)n 

dRd       R  dR    n-rr  AL- 


L 

nrt 


when  R-«-0 


Xi-(-Dn]- 


T^  when  R=R; 


It  is  obvious  that  contributions  to  the  general  solution,  which 
tend  to  infinity  at  (0,0),  must  be  ruled  out  in  these  special  cases. 
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Equations  (98)  and  (9^-),  then,  give 


oo 


2  \       1     d         d       n2^  ^       117x0 

A     A     §  I R  dR       dR       R2$ d     *         $ 

n=1 


n=1  R 

respectively.  For  the  calculation  of  Nu,  Ad  is  chosen  as  «°,   and 


I  2  r22-    did        d      n27^     d  1      ,        n*0 

NuT  >s  -  -   /    ( R 5 r)wAsin 

itA     ^  g  I dR  R  dR       dR         §2     dR  R2     A         § 

n=1 


(108) 
ft-Rj,  RA 


5  2    <~-    1       |        d      nV    .        1 

5,5     §  $2  Z—    H2  dR       dR       2l3$2     *  (-1)n 
n=1  § 


Sirnilar  techniques  can  be  'used  to  solve  equations  (100)   to 
(103)   as  follows i 

2^        n*0 

I 


^c-Y"    sii>3  (110) 


n=1 


where  fis  the  solution  of 


0d2^       dT     n2^  x      n-nr  m     -, 

r2— +  a _.=  .R2-j rs-(-1)nTf  (111) 

dR  dR         <j>  £  L  J 

subject  to  the  boundary  conditions 


T«Tj  when  R=Rn 


■Tj  when  R=R^i# 
Equation  (111)   is  a  Gauchy's  linear  equation.  Its  solution  is 


(112) 


3^ 


H7T  IL7\ 

rmQ^jf    +C2R     *     +P1(R,n)  (113) 

where  P^(R,n)  is  the  particular  integral.  When  tjf  t*  and  Q  are  all 
constants,   the  right-hand  side  of  equation  (111)   becomes 

R2<#r  .    >«]    nTTi 


R2Qjri-(-1)n]       nrtfxL  -(-1)n^ ]r2  n* 

and  P  = — -  -  L-2 — -— U ( — /2) 

1         mxih-n2^/^2)         £(lf-nV/$2)  $ 

R2lnR .  r  7     R2lnR  r  ,  Mx 

C-j   and  C2  can  then  be  easily  determined  by  using  equation  (112), 

When  Rj=0,   C2  must  be  zero  to  avoid  infinite  velocity  at  R=0, 
In  case  ^r%=^=0  and  Q=constant, 

Q$[l-(.1)n]    r  _  fl        -.  jm 

m       TL  '      (D/2)2(2R/D)&     -R2  (_  #>) 

nTK^n2^/^2)  L  J  $ 

Q[l-(-1)n]   ,     D  nTt 

=  — iR2ln~  (—  =2) 

8  2R  § 

For  the  calculation  of  Ku, 

R-+IU     i       2r— 22 — .  R +R.  — - 

T(.XJL    I)=  J\  f(at  R=  JLJL).(.i)   2 

2         2       ^^TT^  2 

n=1,3,b 

is  chosen  as  Tr".  Thus 


si. 
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©O   ,z!~,  *    I        oo 


N%  1=  I   >    —  sin /  \  ^<at  »  -3-A).(-D   2         (11*0 

*  n=1         P  A         V     n=1,3,5 

*U  *  I  8  f"r£  1  n/>^(at  R=  XlWoT'               (115) 

$      n=1  $  /    n=1,3,5 


After  nr  is  obtained,  co  can  be  found  by  the  same  method  from 
equations  (102)  a  -id  (103).   ^e  result  is 

2r~   J        n7# 


co=  —  >     iAjsin— —  (116) 


n=1 

where  co  is  the  solution  of 

dcoidtoln^.ir  „7         jr—  t- 

_ ?  + .^^^PL^jn.  cJ     -c2R     $     -P-        (117) 

dR         R  dR       R       §  Mr  I  J 

subject  to  the  boundary  conditions 

w=o  when  Mh  or  i^  ;  (118) 

i.e., 

nr\  n?r 

a>=C3R      +ClfR  +P2(R,n)  (119) 

where  P2  is  the  particular  integral.   C-  and  C^  may  be  determined 
from  equation  (118).  When  %f  t    and  Q  are  constants, 


p  =  m  _J ±13  .  12 H     5  + 

nix  nix 


fr[l,(,Qn]R' 


nlo1*' 


-g-^X— £--16) 


$' 


{[l-(-1)n]F-t  +(-1)n^}R2         nrr 
-    (  ^   ^     i  f (-—  7*1,    2, 

f 


or  if) 


\  ,    c0  qFi-C-D^r^     ([i-(-i)n]p-^+(-i)nrsl^ 

1r3 £RinR+  — —  +  -* 2 fti- 

2  »+5  3 


8 


n?v 
(—  =1 ) 


(T.      G         Q[l-(-1)n]rf*lnH       {[l-(~1)njF»T^  +  (-1)n^}a2lnB 
12     "  k  96 


(—=2) 


fir 6     °2r-2     ^[1-C-1)nJRlf:LnH       {[l>(>1)n]F-^  +  (-1)ntJy 


20  12 


W 


nrr    . 
(— =if). 

When  Rj=0,  G0  and  C^.  must  be  zero  to  avoid  infinite  temperature 
and  velocity  at  R=0, 
Finally, 

G=  -y~T-[l-(-1)nJ\    ^R&r).  (120) 

n=1  Rj 
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Pipes  of  Other  Cross  Sections 

Method  pX  Thin-Plate  An.aJ.p^As  pointed  out  earlier,  analogy 
to  the  deflection  of  elastically  supported  plates  exists  when  Aj*Q. 
Unfortunately,  this  special  field  of  elasticity  is  no  richer  than 
the  field  of  combined  free-  aid  forced-convection  pipe  flow.  Kow 
solutions  can  not  be  found  by  observing  this  analogy.  However,  a 
special  case  of  the  rectangular  pipe  solution  is  compared  with  its 
counterpart  in  thin-plate  theory  in  the  following  just  to  show  a 
little  detail  of  the  analogy  and  to  check  the  solution  obtained 
by  Fourier  transform* 

Fro::  equation  (30), 

.sin— isin— /    (121) 
IT 


/      /        r      m2      n2 


=1,3,5  n=1,3,5  ran 


-   a2   b2    AJ 


when  FA=0,  Q=constant  and  ^-^"Iky*^"0*   In  other  vords,  equa- 
tion (1?1)  is  the  solution  of  the  equation 
d2    d2 

f(6?  +  of >H+BaAci-% 

J  a^=o  when  F=0,  ?=a,  tf=0  or  ^b 

d2    d2  _ 

l(_  ♦  — )ci=0  when  ^=0,  ^a,  f0  or  y-b 

which  is  exactly  the  same  as  that  governing  the  deflection  of  an 
elastically  supported  rectangular  plate  under  uniform  load,  with 
four  edges  sinply  and  rigidly  supported  (Tlmoshenko,17,  j^7).  A 
careful  comparison  checks  the  correctness  of  equation  (121). 
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Again,  referring  back  to  equations  (21)  and  (23),  equation 
(81»),  with  Q=0,  is  the  solution  of 
d2    d2 

\-0  when  ?=0,  £=a,  ^)=0  or  0=b 

d2    o2 

( — -  +  — rr)^A=0  when  £=0,  pa,  r?=0  or  fi=b 
dp   d^  A       >    b    I      ' 

which  is  the  same  as  that  governing  the  deflection  of  an  elastical- 
ly  supported,  uniformly  loaded,  rectangular  plate  with  edges  simply 
and  rigidly  supported.  The  results  in  the  two  fields  check.  This 
latter  way  of  applying  the  analogy  is  better,  since  PA  is  not  nece- 
ssarily zero,  which  means  a  broader  scope. 

In  cases  where  A=0,  thin-plate  analogy  is  more  fruitful.  For 
erample,  pipes  of  equilateral  triangular,  or  right-angled  isosceles 
triangular  cross  sections  can  be  treated  by  borrowing  known  solutions 
from  thin-plate  theor;,  (Ti.  oshenko,  17,  &57)   as  follows. 

Equilateral  triangle  composed  of  sides 

rm  • ■'#  •/$ 

1  2 
y  =  -(-a-x) 
|J5  3 

1    2 

y  =  -(x-  -a) 
J3    3 

with  tp=F=0  and  Q=constant: 
where  d=2a/3. 
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Right-angled  isosceles  triangle  composed  of  sides 
•  x  =o 
y  =0 
y  +  x  =a 
with  "tw=F=0  and  Q-constantj 

«  ~       -ss— *i  sin *-sin *- 

1  +  J?/2       1+tJ2/2 

n(n2-m2)(n2+n2)2 
r-1,3,5  n=2,4,6 

21-  , 22_u  sin L-sin_JL 

1+^2/2       1+f|2/2. 

n(n2-n2)  (n^n2)2 J 
m=2,4-,6  n=1,3»5 

where  d=2a/(2+//2). 

Temperature  difference  and  useful  parameters  can  then  he 

obtained  as  before. 

It  must  be  noted  that  in  applying  thin-plate  analogy,  genera- 
lity of  the  problem  is  severely  restricted. 

Other  Analogies  p£  Techniques.  For  cases  where  A=0,  the  pro- 
blem can  be  solved  by  using  equations  (30),  (31 )»  (29)  and  (27)« 
Owing  to  the  special  nature  of  equations  (30)  and  (29),  a  number  of 
analogies  such  as  sembrane  analogy,  torsion  analogy,  various  hydro- 
dynamical  analogies  and  heat-conduction  analogy  exist  whenever  the 
condition  suits.  Method  of  complex  variables  is  also  applicable  to 
these  special  cases. 
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OOCBBSfZOHS 

There  are  many  directions  in  which  the  present  investigation 
may  be  extended. 

Firstly,  of  course,  numerical  calculations  including  the  solu- 
tion of  equation  (105)  can  be  performed  and  experimental  work  under- 
taken. 

As  to  the  further  extension  of  the  theoretical  exploration, 
the  author  is  of  the  opinion  that  the  till  for  the  elliptical 
pipes  and  annuli  (at  least  when  Q=0)  can  be  solved  in  terms  of 
Kathieu  functions  by  using  elliptical  coordinate  system.  To  this 
end,  a  treatise  (8)  and  a  paper  (9)  by  LcLachlan  will  be  helpful. 
Relaxation  and  other  finite  difference  methods  may  also  be  applied. 

In  the  present  investigation  wall-temperature  distribution 
is  given.  The  problem  may  be  extended  to  cover  the  cases  where 
heat  flux  at  the  wall  is  prescribed.  Generally,  this  boundary 
condition  is  more  difficult  to  handle.  The  three  known  parameters, 
q,  Q  and  g,  must  fulfill  an  integrated  form  of  energy  equation  in 
this  type  of  proble:  . 

Following  Ostrach,  the  effect  of  friction  heating  and  the  pro- 
blem of  unstable  flow  can  be  investigated,  at  least  for  circular 
annuli.  The  numerical  calculations  involved  will  make  a  conputer 
indispensable. 
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Ananrpow 

Untill  aft     i  final  copy  was  prepared  in  November,  19>8, 

the  author  was  not  axrare  of  the  following  work  si 

Ilan,  L.  S. 

Laminar  heat  transfer  in  rectangular  tubes  with  combined 
free  and  forced-convec+ion.  Anor.  3oc.  Naval  Eng,,  Jour. 
67:  163-167,  Jan.  1955* 

Han  L,  S. 

Laminar  heat  transfer  in  rectangular  channels.  AS1E  paper 
No.  5S-A«12^-.  to  be  ^resented  at  the  annual  meeting  of  the 
ASHE.  Nov.  30-Dec.  5,  1958. 

In  the  first  one,  heat  generation  was  not  considered.  In  the 
second  one  the  effects  of  heat  generation  were  evaluated.  Both 
were  solved  by  using  the  classical  Fourier  method,  i.e.,  Fourier 
series  with  undetermined  coefficients.  The  results  checked  those 
obtained  In  the  present  investigation  by  employing  the  Foiirier 
sine  transfer  . 

The  first  work  mentioned  above  is  known  to  the  author  only 
through  t  9   second  ono. 


the  adthc 
november  28,  1953. 
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Appendix  I 

Nomenclature 

Dimensions i 

m 

Force  (=fliL/t2]). 

w 

Length, 

H 

Mass, 

M 

Thermal  energy  (=[FLJ  =  [KL2/t2]). 

H 

Temerature, 

w 

Time, 

Symbols: 

A 

z-component  of  the  temperature  gradient,  [t/l]» 

a 

bide  length  of  a  rectangular  cross- section,  [l]# 

a 

a/d,  dimensionless. 

B 

Gas  constant,  [FLAT]. 

b 

Side  length  of  a  rectangular  cross  aection,  [l]. 

b 

b/d,  dimensionless. 

bei0(p 

Eei  function  of  the  zeroth  order. 

bei»(y) 

d[bei0(p]/d$-. 

ber0(y) 

Ber  function  of  the  zeroth  order. 

ber0(y) 

d[ber0(r)]/d7. 

C 

A  constant. 

°m 

Constants, 

c 
P 

Specific  heat  capacity  at  constant  pressure,  [q/TM], 

*• 

Specific  heat  capacity  at  constant  volume,  [Q/TM], 

D 

Diameter  of  a  circle  or  sector,  [l]  § 

k6 


d        Hydraulic  diameter  of  the  cross  section,  [l], 

F       N  >  dimensionless  pressure  gradient  parameter. 

/  cpKTr0  y 

/Pr-d^  II 
F.       J»>  — |  dimensionless  pressure  gradient  parameter, 

r  Vaa  ^ 

fz       Body  force  per  unit  mass,  [L/t2]. 

1  Pfz 
G        —  — — g,  dimensionless  mass  through- flow. 

K  k 

g  Mass  through- flow,    [&/t]» 

1     flf 
G»  —  -— g,   dimensionless  mass  through-flow. 

KA     k 

00(p  E0(.ij-). 


Gr       f— i — lyw  ■■■,  Grashof  number. 


GrA         ,  ? •  modified  Grashof  number. 


>2/?f  TP  d3 

grad  Three-dimensional  gradient  vector. 

h  Heat  transfer  coefficient,   q/e° ,    [Q/tL2T], 

i  film 

IQ(^)     Modified  Bessel  function  of  the  first  kind,  zeroth  order. 
Jo^p     Bessel  function  of  the  first  kind,  zeroth  order. 

/3f  d 
K  7Ra(' ),   dimensionless  parameter. 

°p 

k  Thermal  conductivity,    [Q/tL2(T/L)]  . 

KA  "/RaAC         ),   dimension! esc  parameter. 
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w 

Modified  Bessel  function  of  the  second  kind,  zeroth  order. 

kei0(p 

Kei  function  of  the  zeroth  order* 

kei0($-) 

d|]kei0(y)]/d^ 

ker0(p 

Ker  function  of  the  zeroth  order. 

ker^) 

d[ker0(y)]/d5-. 

1 

Half-length  of  the  pipe,  [l] • 

1 

Half-length  of  that  portion  of  the  pipe  vhere  the  flow 

is  fully  developed,  [l]. 

m 

1,  2,  3,  etc. 

N 

Inward  normal  coordinate  of  [* ,  [l]» 

n 

1,  2,  3,  etc. 

Ru 

hd/k,  Husselt  number. 

9 

Pressure,  [F/L2]# 

*« 

Particular  integral  of  a  nonhomogeneous  ordinary  linear 

differential  equation. 

Pr 

<uc_/kf  Prandtl  number. 

Q 

Heat- source  tern,  [Q/tlA]. 

3 

Qd2/kT  ,  dimensionless  heat- source  term. 

4 

Qd/kA,  dimensionless  heat-source  term. 

* 

Heat  flux  at  the  wall,  positive  in  the  direction  of  N, 

[Q/tL2]. 

R 

r/d,  dimensionless. 

r 

Two-dimensional  radial  coordinate,  £l]» 

Ba 

Pr-Gr/7,  Hayleigh  number. 

H*A 

Pr-Gr^/y,  modified  Hayleigh  number. 

S 

Two-dimensional  region  investigated. 

k8 

SA 

RaA^,   dimensionless. 

T 

Temperature,    [TJ. 

T* 

T-Az,    [T]. 

u 

Internal  energy  per  unit  mass,    [q/m]. 

V 

Specific  volume,    [lVm]. 

V 

Axial  velocity,    [L/t]. 

X 

Two-dimensional  rectangular  Cartesian  coordinate,    [lJ« 

*1»   x2 

Two-dimensional  coordinates,   [l]   or  dimensionless,   de- 

pending on  the  specific  coordinate  s/stem  used. 

y 

Two-dimensional  rectangular  Cartesian  coordinate,   [li]. 

*o<p 

Bessel  function  of  the  second  kind  (Weber),   zeroth 

order. 

z 

Cartesian  coordinate  in  the  direction  of  negative  fz,   [lJ. 

bt 

b/a,   dimensionless. 

f3 

Thermal  coefficient  of  volumetric  expansion,    fl/TJ. 

r 

Boundary  of  S. 

To 

A  reference  point  on  f\ 

•y 

Cp/cv,   dimensionless. 

V* 

Laplacian,    [l/L2J. 

T2 

(V2)d2,   dimensionless. 

V* 

VV,    [l/I>]. 

* 

(V  )d%   dimensionless. 

i 

1-^/F,   diiensionless. 

(k 

1-^A7\/FA>   dimensionless. 

s 

A  parametric  variable. 

7 

y/d,   dimensionless. 

$ 

T-TPo,    [T]. 
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•°  A  certain  reference  e,   [t], 

A  D^/Dj9   dimensionless. 

jL(  Dynamic  viscosity!    [M/Lt], 

V  H/d,  dimensionless, 

\  x/d,   dimensionless, 

II  dp/dz  +   pr  fz,   pressure-grdient  parameter,    [f/L^], 

p  Density,  [m/1.3], 

CT  Area  in  the  region  S,    [L2]. 

•"  fr/d2,   dimensionless. 

t  e/Tj.    ,   dimensionless  temperature  difference, 

T^  ©/Ad,   dimensionless  temperature  difference, 

t*  e°/Tp   ,  dimensionless, 

tj[  ©°/Ad,   dimensionless, 

§  Angular  measurement  of  a  ring-  or  full  sector. 

0  Two-dimensional  angular  coordinate,  dimensionless, 

yt  )  A  function  of  one  or  several  variables, 

v 

oo  i     I,       ,   dimensionless  velocity, 

^pKTjr 


Pr 

v 

^  -sBB-asseag  ,   diriensionless  velocity, 

cpKAAd 

Pr 

A  and  Aja  Various  determinants,  dimensionless. 

Subscripts: 

m        1,  2,  3,  etc, 

n        1,  2,  3,  etc. 
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1 

On  the  inner  wall  of  a  circular  annulus  or  ring-sector. 

*1 

n\  -component. 

*2 

xp-component. 

r 

On  r. 

To 

At  r0. 

A 

On  the  outer  wall  of  a  circular  annulus  or  sector. 

* 

On  the  wall  0=§. 

1 

On  the  wall  £=0.^ 

2 

On  the  wall  ^0. 

Applied  to  dimensional  or  dimension- 

3 

On  the  wall  F=a.  \  leas  temperatures  and  temperature 

t 

On  the  wall  ►?=!>• 

differences  only. 

5 

On  the  wall  0=0 J 

» 

"Or";  e.g.,  "I,A,,=  "on  either  the  inner  or  the  outer  wall". 

Superscripts: 

(n) 

i   ,  euc. 

f^f 

Single  finite  FOurier  sine  transform  with  respect  to  F • 

s<- 

Single  finite  Fourier  sine  transform  with  respect  to  n, 

or  0. 

<& 

Double  finite  Fourier  sine  transform. 

\^» 

A  special  finite  Hankel  transform.  (See  Appendix  III). 

Miscellany! 

I 

+,  etc. 

A 

"Plus  (minus)  when  on  the  inner  (outer)  wall",  etc. 

5 

1 

*(-i)n 

i 

"Times  1,  when  on  0=0;  tii.es  (-1)n  ,  when  on  0=$." 
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Appendix  II 

Derivation  of  Equations  (51)  to  (58) 

To  solve  equation  (**1)f  which  is  linear,  the  complementary 
function  or  the  solution  of  the  equation  (Vl)  with  £A=0  is  found 
first  to  be 

C^J0(iisAH)+C^y0(i*sAR)+C^I0(i*sAR)+C4K0(i^sAR) 
(McLachlan,  7,   P.  157,   Formula  7) 
=C^J0(iisAR)+C^Y0(i^sAR)+C^bero(sAR)+iG^beio(sAR)+C4ker0(sAR) 

♦iC4kei0(sAR) 
(McLachlan,  7t  P. 168,  Formula  159;   P. 170,  Formula  201). 
Now,  K0(i-isAR)=ii^[l0(i-isAR)+lY0(iisAR)] 

(LcLachlan,  7,   P. 16^,  Formula  110), 

c«y  (i*sAR)=icii  (i-*sAR) £  K  <r*s.R). 

2oA  20  A  tvO  A 

But  I0(i-isAR)=J0(i*sAR) 

(McLachlan,7,   P.162,  Formula  81), 

1  1  1  2iC» 

C« J  (lis  H)^C«Y/N(itsi.R)=(C»+C'i)J  (its  R) 2.  k  (i-*s.R) 

10  A  20  A  120  A  7ro  A 

=C?  Jo(i£sAR)+CJ|K0(i-*sAR) . 
Again  J0(i^sAR)=J0(i-3/2SAR)ssbero(SAR).lbeio(SAH) 

(   cLachlan,   7,   P.162,  Formula  61 5   P. 168,   Formula  160), 
and  Kq ( i'^s ^R) =ker0 ( s^R) -ikei0 ( s  aR) . 

Therefore,   the  complementary  function  is 

(Cf+CJ)bero(sAfl)+i(CJ-Cf)beio(sAR)+(0i5.+  C5)kero(sAR)+i<C^.C2)keio(sAR) 
=C1  ber0  ( sAR)-»-C2bei0(  sAR)+C3ker0(  sAR)+Cifkei0(  sAR) . 
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When  equation  (50)  holds,  it  is  easily  seen  that 

<^A=C1  ber0(  sAR)+C2bei0(  sAR)+C3ker0(  sAR)+<tycei0(  sAR)+  *A/sA, 

d  r  -,      1     di. 

— HsssJC-1ber,(sAR)+G  bei,(s.R)+C,ker,(sAR)+C,kei«(s,R)  +  mm  —a, 
aa  A     A|_  1        o     A         2       o     A         3       o     A         H-       °     A    J      s^  dR 

A 


r  — w=s     c  Rber*(s.R)+C0Rbeil(sAR)+G-Rker,(s,R)+C,,Rkei,(sAR)l 
ah  A     Ai   1  oA         2         oA         3         oA         *f         oAJ 


1         d    . 

s7       dR  A 
A 


d        d 

—  R  - 
dR       dR 


H — co=s?LcRbei   (sAR)+C0Rber  (s.R)-C.Rkei   (sAR)+C,  Rker  (s.R)l 
dRAALl  o     A         2         oA         3         oAm-         o     A    j 


1  d    d 

+  _ r  _£ 

<J  dR   dR  A 

A 

(KcLachlan,  7,  P. 170,  Formulas  189  and  190;  P. 172,  Formulas  23^ 
and  235), 

1      1   d         d      . 
♦  -_-( h  — )^  . 

s7  R  dR       dR     A 
A 

Boundary  conditions   (h2)   and  C+3)   nov;  beco   e 

•  if 
C1bep0(sARI)+C2beio(sARi)^3kep0(sARi)+alJcelo(sARi)=  -VsA 

.  if 
C1  ber0(  sARA)+C2bei0(  sARA)+C3ker0(sARx)+Clfkei0(sARA)=  -  %/sA 

C1  bei0  ( s  ARj  )-C2beP0  ( sAR!  >+<J3kel0  ( s  Aaz  )-0lfkero  (  s^ ) »  -FA/sA 

C1bei0(sARA)-G2ber0(sARA)^G3kei0(sARA).GlfkeP0(sARA)=  -  *k?k/s\. 

C-j,  C2,  C3  and  Glf  can  then  be  solved  to  be 

C2=^/a 
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G3s^yA 

where  A  and  A^   are  determinants  defined  in  equations  i'^k)   to  (58)« 
Thus,  equation  (51)  is  obtained. 
Furthermore, 

1  d    d 

; h- 

R  dR   dR 


W<-  rr  R  -t^a. 


1   d         d 
and  (—  —  R  — )co    has  already  been  expressed  in  terms  of  bei,   ber, 
R  dR       dR     A 

kei  and  ker  functions;   equation  (52)   follows  directly. 

Finally, 

J* 

GA=ZM     c^RdR 


RI 

.RA 

A 


~  ~AI     L  1  RbSro  (  SAH)+ A2Rbeio ( SAH)+ ^^^o ( 3Aa)+\Hkeio ( sAR)l dR 


RI 


Ari 

which  leads  to  equation  (53)  when  integration  formulas  (KcLachlan, 
7,  P.  170,  Formulas  189  and  190;  P.  172,  Formulas  231*-  and  235)  are 

used. 


^ 


Appendix  III 
Derivation  of  Equation  (60) 
Sneddon  (16)  has  proved  that  if 

f»l  HKR)  [j0(fmR>0o<{'mV-;ro(UB^cio(rBiH)]|M 

JRi 

In  which  5"  &  is  a  root  of  the  equation 

*4fAW&\M*%*\WSA>**  (A1 > 

the  following  inversion  formula  holds  true: 

V*YSm     p  *°     ,Rl?t rljo(^H)Q0(fiaRA)-J0(rmRA)G0($-mR)],(A2) 

L        Jo(^ir^^Jo^mRl) L 

■ 

where  the  summation  extends  over  all  the  positive  roots  %m   of 
equation  (A1), 

Sneddon  (16)  has  also  shown  that 

RdR  da'  j^jyij)   x 


which,  when  applied  twice,  yields 


1  d    d  0    J  (f  B,)  1  d 


d 
1  dR   dR     JQ(r  Rj)  R  dR   dR 


1  d    d 
-( R  — TV 

Rj  R  dR   dR 


Ra 


Jo^imRIJ 


Thus  equation  (U-1)  can  be  transformed  into 
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and 


LvfmRi> 


-<l 


V"* 


K.rVW ., 


Substituting  the  above  expression  for  co^  into  equation  (A2), 
equation  (60)  is  established. 


Appendix  IV 


Derivation  of  Equation  {77) 


It  can  be  shown  (Sneddon,  16)  that  if 


0 


then, 


m=1 


a 


Similarly,  if 


Y(m,n)=\  yCmj^sin-rZiy  , 


then, 


PO         00 


k  r3    r~*     m7t5 

Y*  —  )        )   Ysin~T"sil 


ab 


(A5) 


tf         n7r»7 

—sin— — '  • 
b 

121=1    BP"-1 

Tables  of  y  or  "^  for  simple  functional  forms  of  yean  be  found  in 
Sneddon  (16)  or  Scott  *15).   Repeated  use  of  the  table  gives  corres- 
ponding  values  of  Y.  It  is  only  necessary  here  to  mention  that,   for 
a  constant  C,   the  following  finite  Fourier  sine  transforms  are  list- 
ed also  in  Scott  (15)   and  Sneddon  (16): 


O 


cm 


[l-(-1>B] 


C=  —c[l.(-l)n] 


(A6) 
(A6«) 
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and 


*       at) 


-C  [l-<-1  )m]  [1-C-1  )a]. 

ronTT 


It  is  not  difficult  to  show  (Sneddon,   16)   that 


ddY     m*, 


»v 


And  it  is  still  easier  to  see  that 


a 


62y  d2y 

d^2   d^2 


Therefore. 


-V 


d2T 


which,  when  applied  twice,  yields 


d2  d2 


(A7) 


(A8) 


(A9) 


(A10) 


(— *  +  — r)2^  H  (—5-  +  — o> Y  -(-1  )m|  <— o  +  — s)Y 

dp       d^2  a    Udj2      o>]2     J(0,n)  »■  dj2      o^j2     J 


(*•')). 


-  -=g-[Y<Of«0).(.1)nTKaf7)]+  -^pY 
a  a 

fe^S   d2Y       BITf  d2    r 

-2—-. + YC0f0)-(-1)oytg$n) 

§2     drj2       a     dr|2L       '/  '/. 

In  case  Y(0,^)="Vta,r?)=0,  equation  (A11)  becomes 


d^y 


**' 


(A11) 


d2        d2     0        mTrrr  d2         d2 

(—2 +  t?)  2^=  H  (t  +  — ?>n 

dj2      d»p  a  U  dp      d^j2    -KO,^) 


tfl 


.(-1)-  (— sr+  — 5>)^ 

L  o£2      d^2   J(aVjO 
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nM*^     m2^  d2Y      d^y 

+  — r-Y-2 +  — »  . 

a4"  ap     do2       d<|f 

Analogous  to  equation  (A8),   there  is  the  equation 
d2f      nTT^  ~         _n      n2*2  ^ 

Xj|  =  -\f(C)o)-(-i)^,b)].irr 

which  yields,  when  applied  twice, 


(A12) 


(A13) 


d^      n^rr  d2T 


d7      b 


(_)  -(-1)nK)         _    J 

L  d72  (w=0)  d^  (k7=b)  J 


Now,  if  Y(J,0)=X^,d)=0, 


4*y 

^  (vj=0)  d»j2  (*y=b) 

n3fi3r^  ^^       _-]       n^7r*^ 

[TKmf0)-(-1Jl'Y<m>b)J+  -^V 


a 

f(„o)=r^(^7)gi^](?=o) 


(A1VJ 


i 


ffi?Xf 

Y<£tO)siir^dj: 


CA151) 


~°1 


and,   therefore, 


(A15) 


d*> 


df  *  "     b2    Y 
dH'      n^rr  d2T 


"T- 


<—*> 


d2f 


*f> 


-C-1)a(— W 


1     n  7r    ^ 


drj*       b    U^j2   ("|=0)  dn2   (^=b)-l        b 


(A16) 


(A17) 
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Furthermore,  if  "f(0,n)«y(af  »7)=0,  equation  (A10)  becomes 


r.V„     d2f 


-y+ 


.2' 


kd^2      d^2'  a"  do" 

which  gives,  vhen^C  FtO)=fy<J,b)=0     noting  equations   (A15)   and  (A15')> 


d2f 


2     J~ 


(0=0)    LV2      d  : 


f  j  (? 


(1=0) 


d"Y 
(-4) 


a2      a2 


d*p  (K|=b)    L  d^       b*f    J(v?  = 
Hence,    equation  (A17)   becomes 


d*V      nTrfr  ^2        *2 


— r 
drj* 


b  ILa^2     a*?2   J(w=o)  U52     d72-i(^=b)J 


n^TTsst 

b* 


(A18) 


Using  equations  (A16)   and  (A18),   equation  (A1?)  is  transformed  into 


a2      a2   0     urn  rr  a2      a2     7  ra2      a^ 

df2       d*P     r     a   Ud^2       df     %»0>  Ld^2       ar|2    J(£=a) 


a2 


ra2       a2 


_)y]        -(-i)n(-r+  -*>* 
L  dp     a*?2    J(n=o)  L  dj2     a^j2    J(y=b) 


,    m^       n^      ^ 
a         b 


(A19) 


Now,  coA(0,^?)=^(a,y)=^(c,0)=WA(F,b)=0.   So  equation  (63)   can  be 
transformed  by  using  equations  (A19),    (70),    (A6),    (A6')   and  (A7) 
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into 


2         2 
mw       n^ 


\    1a.       V  r-  ~rr-  _      Dlb  8J1 


a1 


an         b 


a    -  --     b 

As 

Substituting  wA  obtained  from  equation  (A20)  into  equation  (A5), 

equation  (77)   is  established. 

Similar  method  can  be  used  to  solve  equations  (73)  to  (76) 
yielding  equations  (flfc)  and  (  V). 


Appendix  V 


Derivation  of  Equation  (10^) 
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Set 


As 


o 


rojb 
0)sin — -d0 


then,   si  ilar  to  Appendix  IV, 
R  dR       dR       R^  60^         6R^ 


R  dR       R  00' 


d2^      1  dy     1     n2rr2^    1 


+  — 


dR£ 


R  dR 


R2     $2         *■  $ 


•A  -5— [■»KRf0)-(-1)nV<Ht#)J 


and 


iV 


2  dJ?       1 


( »-*«<f  — *)^  -tt  ♦ f  *  -r(1+ 

R  dR       dR       r  d0^  dR         R  dR-3       R*  $ 

1  2n2?<2  d^      1     n2^2  nV         a, 


2n2^2  d2r 


r3 


^     dR 


2     ^«        «*+     *2 


(- 


¥     9 

1     n?rr1   d         d         1      d2 


-W 


R2  §    i-  R  dR       dR       R2  602 

2 


(R,0) 


1  d 

-(-1)n(-- 

R  dR       dR       R  d0^ 


6        1  d 


(R,$) 


+  (• 


1   d 
+  — 


dR' 


1     n27^  (1     n-rry 

■>  ~*  -~fytR,0)-(-1  )n^R,o)     . 


R  dR        R«*     ^      IJP  $ 


Therefore,  using  equations  (96)   and  (97), 
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1  d 


1     o: 


d  co. 


d3<& 


2nV1d% 


R  dR       dR       R2  d02       A     dR*         8  dRJ         t*  |" 

> 


$z     dR' 


1  2n2-rr2  dco.       1     nV  nV 


I 

1       TWr  fr 


$2    'dR  R*     §2     "    $ 


r(H-1)n]FA"V(-1)ntA$} 


R*  § 

and  equation  (95)  can  be  transformed  Into 

d**^       2  d3w        1           2n2^  d2&       1           2nV  dco 
Jk  + Jk -d+  .—.) — «4  ♦  -—.(1+  z~) — • 


dB 


R  dR3         fl2 


dR' 


*-x 


H"L  f       $■ 


n3 


1     mrr,- 


$2     dR 


( — r — l+)+Ra. lot  *JL-  -—  — 


A    A     A 


E*  $ 


'-(-1)n]FA-V(-1)ntA*}- 


This  variable-coefficient,  linear,   non- homogeneous,   ordinary  differ- 
ential  equation  can  be  solved,    subject  to  the  boundary  conditions 
kjy»0  when  R^Rj  or  R 

^  d2&       1   d<& 


1  d         d         1     d 

R  dR      dR       «T  d0^     A     dR^ 


+   _  — • 
R  dR 


'  ZFA[1-(--nn]-1ii  when 


n7r 

| 

n^ 


*A[1-C-Dn]-^X  when  R«RA, 


Hence,  equation  (ICh)   is  established. 


A  THEORETICAL  INVESTIGATION  OF  COMBINED  FREE-  AND  FORCED- 
CONVECTION  HEAT-GENERATING  LAMINAR  FLOW  INSIDE 
VERTICAL  PIPES  WITH  PRESCRIBED 
WALL  TEMPERATURES 


by 


PAU-CHANG  LU 
B.  S.,  National  Taiwan  University,  China,  1953 


AN  ABSTRACT  OF 
A  THESIS 


submitted  in  partial  fulfillment  of  the 


requirements  for  the  degree 


MASTER  OF  SCIENCE 


Department  of  Mechanical  Engineering 


KANSAS  STATE  COLLEGE 
OF  AGRICULTURE  AND  APPLIED  SCIENCE 


1959 


Recently  the  free- convection  effect  in  vertical  pipe  flow 
has  become  quite  significant  in  the  heat- exchanger  problem  of 
reactor  engineering  and  the  blade-cooling  problem  in  gas  turbine 
engineering.  The  modifier  "vertical"  refers  to  the  general  assump- 
tion that  the  longitudinal  dimension  of  the  pipe  is  parallel  to 
the  direction  of  the  generating  bod/  force  such  as  gravity  or 
centrifugal  force. 

It  is  the  purpose  of  the  present  investigation  to  find  the 
velocity  and  temperature  distributions  inside  vertical  pipes  of 
various  cross  sections,  through  which  heat-generating  fluids  are 
flowing  laminarly.  Analytical  means  were  employed  throughout. 
The  governing  equations  of  the  problem  were  established  first  in 
dimensionless  forms  subject  to  nine  general  assumptions.  The  most 
important  assumptions  are: 

1.  The  fluid  is  semi- incompressible. 

2.  Energy  dissipation  is  negligible. 

3.  The  axial  component  of  the  temperature  gradient  is  constant. 
With  the  prescribed  wall- temperature  distribution  and  the  known 
mass  through- flow,  the  equations  were  solved  for  a  number  of  spe- 
cific cros-  sections.  Nusselt  numbers  at  the  wall  were  also  formu- 
lated. 

Plow  inside  circular  annuli  was  described  in  tarms  of  ber, 
bei,  ker  and  kei  functions.  The  finite  Fourier  sine  transform  was 
employed  in  obtaining  solutions  for  rectangular  pipes  and  pipes  of 
ring- sector  or  full  sector  cross  sections.  A  special  finite  Hankel 
transform  was  used  to  find,  formally,  the  more  general  solutions 


for  circular  annul!.  These  specific  solutions  are  all  original  con- 
tributions of  the  present  investigation. 

Also,  an  analogy  to  the  deflection  of  elastically  supported 
thin  plates  was  observed  for  the  first  time  in  some  special  cases. 

The  possibility  of  constructing  solutions  for  elliptical  pipes 
and  annul!  in  terms  of  Mathieu  functions  was  pointed  out  together 
with  some  other  possible  directions  In  which  the  present  investi- 
gation may  be  extended. 

Most  of  the  mathematical  methods  used  were  employed  as  formal 
tools  only.  Unyielding  rigour  was  not  attempted.  The  justification 
of  the  results  depends  generally  on  the  successes  of  the  same  me- 
thods in  other  fields  of  application. 

Numerical  calculations  were  not  taken  up.  No  verifying  experi- 
ments were  conducted. 


